Krishna et al. Theoretical Biology and Medical Modelling 201 3, 1 0:54 
http://www.tbionned.conn/content/1 0/1 /54 



THEORETICAL BIOLOGY AND 
MEDICAL MODELLING 



RESEARCH 



Open Access 



Rate-dependent Ca^"^ signalling underlying 
the force-frequency response in rat ventricular 
myocytes: a coupled electromechanical 
modeling study 

Abhilash Krishna\ Miguel Valderrabano^, Philip T Palade^ and John WCIark JrJ" 



^Correspondence: jwc(5)rice.edu 
^ Department of Electrical and 
Computer Engineering, Rice 
University, Houston, Texas, USA 
Full list of author information is 
available at the end of the article 



Abstract 

Background: Rate-dependent effects on the Ca^+ sub-system in a rat ventricular 
myocyte are investigated. Here, we employ a deterministic mathematical model 
describing various Ca^+ signalling pathways under voltage clamp (VC) conditions, to 
better understand the important role of calmodulin (CaM) in modulating the key 
control variables Ca^+/calmodulin-dependent protein kinase-ll (CaMKII), calcineurin 
(CaN), and cyclic adenosine monophosphate (cAMP) as they affect various intracellular 
targets. In particular, we study the frequency dependence of the peak force generated 
by the myofilaments, the force-frequency response (FFR). 

Methods: Our cell model incorporates frequency-dependent CaM-mediated spatially 
heterogenous interaction of CaMKII and CaN with their principal targets 
(dihydropyridine (DHPR) and ryanodine (RyR) receptors and the SERCA pump). It also 
accounts for the rate-dependent effects of phospholamban (PLB) on the SERCA pump; 
the rate-dependent role of cAMP in up-regulation of the L-type Ca^+ channel (/ca,/.); 
and the enhancement in SERCA pump activity via phosphorylation of PLB. 

Results: Our model reproduces positive peak FFR observed in rat ventricular myocytes 
during voltage-clamp studies both in the presence/absence of cAMP mediated 
yS-adrenergic stimulation. This study provides quantitative insight into the 
rate-dependence of Ca^+-induced Ca^+-release (CICR) by investigating the 
frequency-dependence of the trigger current {Ico,l) and RyR-release. It also highlights 
the relative role of the sodium-calcium exchanger (NCX) and the SERCA pump at 
higher frequencies, as well as the rate-dependence of sarcoplasmic reticulum (SR) Ca^+ 
content. A rigorous Ca^+ balance imposed on our investigation of these Ca^+ 
signalling pathways clarifies their individual roles. Here, we present a coupled 
electromechanical study emphasizing the rate-dependence of isometric force 
developed and also investigate the temperature-dependence of FFR. 

Conclusions: Our model provides mechanistic biophysically based explanations for 
the rate-dependence of CICR, generating useful and testable hypotheses. Although rat 
ventricular myocytes exhibit a positive peak FFR in the presence/absence of 
beta-adrenergic stimulation, they show a characteristic increase in the positive slope in 
FFR due to the presence of Norepinephrine or Isoproterenol. Our study identifies 
cAMP-mediated stimulation, and rate-dependent CaMKII-mediated up-regulation of 
lca,L as the key mechanisms underlying the aforementioned positive FFR. 
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Background 

Cardiac muscle contraction is a result of a transient increase in intracellular Ca^^ 
concentration [C/^-^+l^yo. Sarcolemmal (SL) membrane depolarization triggers Ca?'^ 
influx via the dihydropyridine (DHP) -sensitive L-type Ca?'^ channels. Following diffu- 
sion across a small sub-membrane dyadic space, this influx activates ryanodine receptors 
(RyRs) controlling ryanodine-sensitive Ca}^ release channels in the junctional portion 
of the sarcoplasmic reticulum (jSR). Fabiato and Fabiato [1] named the process calcium- 
induced calcium release (CICR). Ca^^ subsequently diffuses from the dyadic space into 
the myoplasm. Ultimately, intracellular Ca?'^ concentration [Ca^~^]myo is returned to rest- 
ing levels by combination of: (a) Ca^^ buffering in the dyadic space and myoplasm; (b) 
sequestration of Ca^~^ by sarcoplasmic/endoplasmic reticulum C<3^^+-ATPase (SERCA)- 
type calcium pumps lining the longitudinal portion of the sarcoplasmic reticulum (LSR); 
and (c) Ca^'^ extrusion from the myoplasm hy Na^ I Ca^^ exchangers and C<3^'^'^-ATPase 
pumps on the sarcolemmal membrane. 

Ca^^ is an extremely important and highly versatile second messenger in cardiac cells, 
which plays a crucial role not only in excitation-contraction (E-C) coupling but also in 
excitation-transcription coupling [2]. Various inter-connected Ca?'^ signalling pathways 
help preserve the integrity of the cellular Ca^^ system despite any changes in pacing fre- 
quency. Specifically, Ca^^ triggers the CaM-mediated rate-dependent effects of CaMKII 
and CaN on the characteristics of the apposed dihydropyridine (DHP) and ryanodine- 
sensitive Ca}^ channels in the dyad, whose interaction forms the basis for CICR. The 
proteins CaMKII and CaN not only influence the release mechanism but also affect 
the SERCA pump either directly or indirectly via phospholamban (PLB), thus modulat- 
ing the Ca?^ uptake process [3-6]. It is also known that ^-adrenergic stimulation of the 
cardiac cell, mediated by the second messenger cAMP, modulates the frequency depen- 
dence of the peak force generated by the myofilaments, the force-frequency response 
(FFR). These protein-mediated and second messenger pathways help maintain Ca^^ 
homeostasis over a wide range of stimulation frequencies. 

Cardiac contractile function is closely coupled with heart rate (Bowditch effect [7]). 
Although positive [8], almost flat [9] and negative [10,11] peak FFR have been reported 
in the literature, it is clear from in-vitro studies involving stimulation in the physiological 
range of frequencies [12,13] that rat ventricle exhibits a positive peak FFR. The issue of 
force-frequency relationship requires broadened investigation into the various underlying 
cellular and molecular mechanisms. We propose that mathematical modeling would be a 
useful tool in helping to sort out this complex issue. 

Methods 

All simulations and analysis were performed on a 2.8GHz Intel® Core™2 Duo 
CPU-based computer using Microsoft Windows XP operating system. The sarcolemmal 
membrane charge balance equations, the Ca^^ material balance equations in the 
myoplasm and SR, and the force balance equations describing the model for myofil- 
ament contraction constitute a set of 93 ordinary differential equations (ODEs). A 
fixed-step Merson-modified Runge-Kutta 4th-order numerical integration scheme [14] 
was used to solve this set of Ist-order differential equations (ODE) describing the dynamic 
model. The free Ca^^ concentration in the dyad is governed by the time courses of the 
Ca}^ fluxes through Ca^^ transport systems, as well as by the time course of Ca^'^ 
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binding to Ca?'^ buffers present in the junction [15]. Description of the spatio-temporal 
dynamics of calcium transients in the dyad triggered by Ca^^ stimulus (basis of CICR) 
requires calculation of the partial differential equations (PDE) of the whole reaction- 
diffusion system. Formation and dissipation of Ca^^ gradients around an open channel 
(DHP-sensitive and Ry-sensitive channels in the dyad) is assumed instantaneous as was 
validated for microsecond timescale and nanoscopic space by Naraghi and Neher [16]. 
Local Ca^^ concentration in the vicinity of open channels (located on opposing bound- 
aries of the dyadic space) was calculated as the steady state gradient around a point 
source [17]. The Ca}^ concentration increments from individual channels at each point 
in space were assumed to be additive [16,18]. The software kernel follows the changes in 
the state of trigger and release channels together with variables like membrane voltage 
and spatial Ca^^ concentration to calculate the instantaneous rate constants and estimate 
the duration of transient events. Crank [19] discusses diffusion problems in a two-phase 
heterogeneous medium and shows that diffusion through a system of barriers (RyR feet 
structures in the dyadic cleft space) can be approximated by diffusion in the same region 
without barriers but with a reduced effective diffusion coefficient. We hence take this 
approach in modeling the CcP'^ diffusion by solving the 2-D Laplacian equation (Krishna 
et al. [15], Appendix A3, Eq. 140) in the DCU without explicitly accounting for local 
potential fields. More specifically, an explicit finite difference scheme was used to solve 
these Laplacian equations describing C<3^^+ -diffusion in the dyadic space analogous to the 
method detailed in Smith et al. [20]. Specifically, a radial symmetry is employed in solv- 
ing the PDE in the dyadic volume allowing the solution to be computed in a rectangular 
cross-section discretized into a 20 by 20 cartesian grid. The spatial step size used in the r 
and z-direction (Figure IB, Krishna et al [15]) was 10 nm and 0.76 nm respectively (Table 
two, Krishna et al. [15]). The 20x20 grid size was used to obtain a stable numerical solu- 
tion using the explicit finite difference scheme employed to solve the PDE. Obtaining an 
accurate description for C^^^"^"*" -diffusion in the dyadic space is vital to ensure adequate 
time delays associated with RyR release (z-direction) and C<2^+ -diffusion into the cytosol 
(r-direction) which controls the rate of SR C<3^^+ -uptake via the SERCA pump. Both 
of these delays are important in ensuring robust luminal sensor mediated RyR refrac- 
tory characteristics (described in Krishna et al. [15]). We use the method of lines 
(discretization in space) to solve the PDE. The full set of ODEs and finite difference 
equations are solved simultaneously to obtain the complete solution. Execution of a sin- 
gle cycle which translates to 200 ms at 5 Hz took 21 seconds with a time step of 1/xs. A 
non-linear leastsquares method [21] was used for parameter estimation and data fitting. 
Results were visualized using Matlab (Mathworks, Natick, MA) and Origin (OriginLab 
Corp., Northampton, MA). 

Model development 

Our objective was to develop a model of the rat ventricular cell under voltage clamp con- 
ditions, which includes the description of various Ca?'^ signalling pathways in the dyadic 
space, the myoplasmic medium and the sarcoplasmic reticulum. Therefore, we start with 
a broad discussion of the key proteins involved in the Ca?^ signalling mechanism and 
continue with a progressively more detailed description of their influence on respective 
targets. It is important to note that all Ca?'^ concentrations discussed in the model pertain 
to unbound Ca^^ unless specified. A detailed description of the membrane classification. 
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(B ) Probabilistic model of CaMKII subunit switching 
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Reaction map for reversible binding of CaM, Ca2CaM and Ca4CaM to CaN 
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Figure 1 CaM-CaMKII-CaN pathways. Reaction Maps. (A) Cooperative binding of 2 C6/^+ to CaM 

sequentially to the C and then N-terminal hands along with binding of CaM buffers; (B) Probabilistic model 
of CaMKII subunit switching; (C) Reaction map for reversible binding of CaM, Ca2CaM and Ca4CaM to CaN. 



channel and exchanger distribution as well as the various fluid compartments involved is 
given in Krishna et al [15]. 

Calmodulin mediates the regulation of a variety of C<3^^+ -dependent signalling path- 
ways in the heart involving CaMKII and CaN [2,22]. These protein-mediated interactions 
form the basis for a robust mechanism that enables a cells response to increased heart 
rate. CaMKII is reported to be responsive when targeted to Ca^^ release sites such as the 
dyadic cleft, and CaN is responsive to gradual changes in the lower-amplitude myoplas- 
mic Ca^^ signals [2] . This heterogenous response is a result of the different affinities of 
CaMKII and CaN for CaM [22] and the non-uniform distribution of these proteins in the 
cell. Recent study [23] also attributes a key role in frequency-dependent acceleration of 
relaxation to activated-CaMKII in the cytosol. The model of the CaM-dependent Ca^^ 
signalling process (Figure 1), which includes a reaction map for cooperative binding of 
Ca^^ to CaM, the scheme for CaM buffering, probabilistic model of CaMKII subunit 
switching and the reaction map for reversible binding of CaM, Ca2CaM and Ca4CaM to 
CaN, is adopted from Saucerman et al. [2] . However, to reproduce relative local CaMKII 
and CaN activity, modifications were made to the rate constants for CaM buffering in 
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the dyad (Table five, Krishna et al [15]). Specifically, a limitation of the Saucerman and 
Bers model is that it is based on little available information regarding the operation of 
CaM buffers [24] at locations where CaM encounters very high Ca^^ concentrations 
(compared with the myoplasm). 

We incorporate the effects of -adrenergic stimulation via cAMP-dependent mod- 
ulation based on a model (Figure 2) derived from Demir et al. [25]. Stimulation of 
^-adrenoceptors by Isoproterenol (Iso) results in the activation of a G protein (g^) that 
stimulates Adenylate cyclase (ADC) and enhances the production of cAMP. Subsequently, 
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Figure 2 cAMP and cGMP pathways. Reaction pathways triggered by sympatlietic and parasympatlietic 
neural stimuli invoking cAMP [25] and cGMP-mediated [26] modulation of various sub-cellular targets 
respectively. ACh-mediated effects of junctional receptor on A/ci+ current and background current and 
extrajunctional muscarinic Mi/^ACh receptor on Ik^ach (direct muscarinic pathway) and cAMP-mediated 
effects of ;0-adrenoreceptor (adrenergic pathway) and cholinergic M2/ADC receptor (indirect muscarinic 
pathway) on L-type Ca^+ (lca,/_)/ ik), hyperpolarization-activated (If), and No'^/K'^ (l/vo/()currents are 
shownf. Although cAMP-dependent modulation of \ca,L and lup are known to be PKA-mediated, we make 
use of a lumped cAMP term to model this influence. A similar lumped cAMP term is employed to modulate 
phosphorylation of PLB which in turn affects the SERCA pump. ACh triggered NO-mediated synthesis of 
cGMP by sGC is modeled employing a 3-state model for sGC transition based on Yang et al. [26]. cGMP is 
involved in suppressing cAMP activity via PDE. cGMP also enhances Irmca, Co^^ activated /(+ channel {Ikco) 
and suppresses Icq,l, Icyt.serca via PKG (not explicitly modeled, but lumped into a cGMP term). fThe rat 
ventricular cell model used in this study is however limited to Cci^+ related channel, exchanger and pumps 
ilca,L, iNaCa, IpMCA and Icyt.serca), while lacking exclusive or /(+ related channels and transporters (as shown 
in Figure 2, Krishna et al. [1 5]). The part of the model describing cAMP-mediated pathway used in our study is 
highlighted (blue). 



Krishna etal. Theoretical Biology and Medical Modelling 2013, 10:54 
http://www.tbiomed.eom/content/1 0/1 /54 



Page 6 of 26 



cAMP may directly or indirectly activate various intracellular targets including ion chan- 
nels and exchangers. The indirect modulation involves activation of cAMP-dependent 
Protein kinase A (PKA) before modulation of the channel protein. The reaction kinet- 
ics for the cGMP-mediated pathway (Figure 2) involving acetylcholine (ACh), nitric 
oxide (NO) and soluble guanylate cyclase (sGC) are adopted from Yang et al. [26]. 
Although it is well known that cGMP modulates its targets via protein kinase G (PKG) 
or phosphodiesterase (PDE), we have refrained from modeling these protein interactions. 
Given that NO synthase inhibition and/or NO donor had little or only marginal effects on 
the FFR in rat ventricular myocardium [27], the cGMP level is kept constant in this study. 

DHP-sensitive Ca^""" channel 

Upon cell depolarization, the L-Type Ca^^ channel {Ica,L) brings trigger Ca^^ into the 
dyad that facilitates Ca?''^ release from the apposed Ry-sensitive Ca}^ channel asso- 
ciated with the membrane of the junctional sarcoplasmic reticulum (jSR). Feedback 
controlled interaction between these apposed channels is critical in CICR as well as in 
the maintenance of the overall cellular Ca^^ balance. Besides membrane voltage, gating 
of the Ica,L channel is also influenced by two prominent C<3^'^"'" -mediated effects, namely 
0/3^-^+ -dependent inactivation (GDI) and 0/3^-^+ -dependent facilitation (CDF). 

In GDI, Ca^^ that enters the dyad either via the Ica,L channel or through RyR release 
from the jSR binds to the protein Calmodulin (CaM) which is tethered to the C- terminus 
of the Ica,L channel [28], modulating the interaction of CaM with the Ca^^ channel. 
Leucine- Alanine (LA) and Isoleucine-Glutamine (IQ) are 2 adjacent motifs in the Ca?'^ 
sensing domain of the C-terminus of the Ica,L channel. A C/s^"^"^ -dependent switch of CaM 
from LA to IQ motif removes CaM from the inner mouth of the channel pore, thus caus- 
ing an enhancement in inactivation by facilitating the constriction of the pore [29,30]. 
GDI is a critical negative feedback mechanism which causes decreased Ca?'^ entry via 
J^Ca,L when the SR load is high with an accompanying large myoplasmic Ca^^ transient, 
and it results in increased Ca^'^ entry via Ica,L when [Ca?'^]niyo is small due to a low SR 
load. We utilize a 2-state model, as shown in Figure 3A of our previous study Krishna 
et al. [15], to simulate GDI. 

In contrast, C/^"^"*" -dependent facilitation (CDF) is a C^^^-^+ZCaM-mediated enhancement 
in Ica,L via activation of CaMKII [31], which has been shown to tether to the Ica,L channel 
[32], functioning as a Ca^^ signalling sensor for facilitation. Activated Calcineurin (CaN) 
is also observed to facilitate Ca^^ entry via the lca,L channel [33]. CDF is implicated 
in causing a gradual increase in amplitude and an accompanying decrease in inactiva- 
tion over consecutive pulses after a resting interval [34]. Enhancement of Ica,L caused 
by activated-CaMKII and CaN plays a key role in negating the effects of incomplete 
ica,L channel recovery at faster heart rates, thus helping to improve cardiac performance 
during exercise. Although CDF and GDI of Ica,L coexist, GDI responds much faster 
(within the same beat) than CDF (over several beats). Our model incorporates CDF by 
allowing the rate constants in the 6-state Markovian model of the Ica,L channel {k^^ , 
and k^^^ in Figure 3A, Krishna et al. [15]) to be a function of the available active 
CaMKII and CaN. 

Rate-dependent increases in lca,L can also be caused by frequency-dependent increase 
in ^-adrenergic stimulation increasing the level of available cAMP (Table 1), which 
in turn causes an enhancement in lca,L channel current via protein kinase- A (PKA). 
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Figure 3 Iryr- Frequency-dependent characteristics of tine SR Cci^+ release cliannel Iryr. Voltage clamp 
protocol used is a 50 ms step pulse to 10 mVfrom a holding potential of -40 mV. (A) RyR open probability vs 
[Co^'^]ryr for increasing frequency of stimulation (0.5 Hz to 1 0.0 Hz). The inset shows the 
frequency-dependent modulation of peak RyR open probability. (B) Frequency dependence of average 
charge transported by NCX and average [A/a+]myo in the myoplasm. (C) Rate dependence of pre-release 
maximum and post-release minimum SR Ca^+ level. (D) Traces for RyR release for the corresponding 
frequencies of stimulation in panel A. (E) Dependence of average level of activated-CaMKII and peak Iryr 
current on frequency of stimulation. (F) Frequency-dependent changes in average Iryr activation and 
inactivation rates. The legend for traces in panel A is shared by panel D. 



Although PKA is involved in the indirect regulation of the Ica,L channel, its effect is con- 
sidered lumped into the conductance term in the ionic current description (Appendix 
A3, Equations 1-2). The effect of ^-adrenergic stimulation via cAMP, particularly its 
dose-dependent influence on L-Type Ca^^ channels both in terms of modifying the 
single channel behavior such as Ca^^ ion permeability as well as overall channel recruit- 
ment characteristics, is not clearly understood. While cAMP has been shown to increase 



Table 1 Frequency dependence of intracellular cAMP concentration 



Stimulation frequency (Hz) 


Intracellular cAMP concentration {fiM) 


0.50-1.67 


2.67^ 


2.50 


3.00 


2.86 


3.19 


3.33 


3.45 


4.00 


3.85 


5.00 


4.67 


5.71 


5.58 


6.67 


18.06 


8.00 


24.58 


10.00 


27.03 


12.00 


29.17 
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various stimulation frequencies used in the model (based on Equation 2, Demir et al. [25]). 
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channel open probability [35,36], increased levels of cAMP also result in increased phos- 
phorylation of L-type Ca^^ channels, causing an increased permeability to Ca^^ ions 
[36-38]. 

Ry-sensitive Ca^"*" release channel 

Ry-sensitive Ca?'^ channels on the jSR membrane respond to the trigger Ca^^ entering 
the dyadic space via the Ica,L channel on the plasma membrane. A larger Ca?-^ release 
from the jSR follows, forming the basis for Ca}^ induced Ca^^ release (CICR) and sub- 
sequent contraction of the ventricular cell. CaM that is tethered to the Ry-sensitive Ca^^ 
channel [39] facilitates frequency-dependent CaMKII and CaN-assisted modulation of 
that channel Although CaMKII is known to bind to the RyR [40-42], the effect of this 
association has not yet been resolved. In lipid bilayer studies, CaMKII has been shown 
to increase [40,41,43] or decrease [44] RyR open probability. In studies on rat ventric- 
ular myocytes, it has been shown that endogenous CaMKII has an activating effect on 
the RyR Ca^^ release channel [45-48]. However, a contrasting study shows that consti- 
tutively active CaMKII depresses RyR release [49]. Thus, the functional consequence of 
phosphorylation of RyR by CaMKII remains controversial. Since the bulk of the pub- 
lished literature on this topic points toward an activating effect of CaMKII on RyR, this 
concept is adopted in our model by making the rate constants in our 4-state Markovian 
model of the RyR channel {l<i2, 1<^^ , l//^ and Figure 3B, Krishna et al [15]) func- 

tions of the available active CaMKII. Although CaN is reported to regulate ryanodine 
receptor Ca^'^ release channels in rat heart [50], we have refrained from modeling its 
influence on the RyR channel because CaN is known to be constitutively active in the dyad 
[2] exhibiting only minor frequency-dependent modulation in its level, hence making its 
rate-dependent regulatory role insignificant. 

SERCA pump 

In rat ventricular myocytes, 85-90% [51,52] of the systolic increase in Ca^^ in the 
myoplasm is recovered back into the SR stores via the sarcoplasmic reticulum Ca^^- 
ATPase (SERCA) pump. Frequency-dependent CaMKII activity is known to cause an 
acceleration of relaxation [3,23]. CaMKII affects the SERCA pump via direct phosphory- 
lation, assisting in enhancement of SR Ca^'^ transport by increasing the pumping rate [4]. 
This feature is incorporated in our model by letting the rate constants for Ca^^ binding 
to/release from the SERCA pump depend on the available active CaMKII. The SERCA 
pump is also indirectly affected by CaMKII via phosphorylation of unphosphorylated 
phospholamban (PLB), relieving the inhibition caused by PLB on the SERCA pump and 
thereby increasing the sensitivity of the pump for Ca^~^ uptake. This indirect effect is 
modeled by having the rate constant for phosphorylation of PLB be a function of active 
CaMKII in the myoplasm. These two effects cause enhancement in SR Ca^~^ uptake in 
an activity- dependent fashion. Expression of CaN, a protein that is a phosphatase, para- 
doxically has been reported to cause an increased level of PLB phosphorylation via an 
unknown indirect mechanism in transgenic mice. In CaN knock-out mice, decreasing 
phosphorylation of PLB allowed an increased level of inhibition of the SERCA pump, 
which resulted in poor muscle contraction and relaxation [5,6,53]. However, it is unclear 
if this behavior is a result of other compensatory mechanisms such as decreased CaMKII 
expression or enhanced PLB to SERCA ratio. On the contrary, CaN has been reported to 
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inhibit SERCA activity in isolated non-failing human myocardium [54] in-vitro. Hence, 
we model the role of CaN in rate-dependent inhibition of the SERCA pump via PLB 
dephosphorylation by allowing the rate constant for phosphorylation of PLB to be depen- 
dent on available active CaN in the myoplasm. With increasing heart rates, -adrenergic 
stimulation results in increasing levels of cAMP in vivo (Table 1 shows values for maximal 
-adrenergic stimulation), which in turn phosphorylates PLB via PKA, causing enhanced 
uptake by the SERCA pump (Appendix A3, Equations 5-6). Together, activity- dependent 
recruitment of these CaMKII and cAMP-mediated effects at high frequencies counter the 
effect of CaN as well as decreasing cardiac cycle duration on SR refilling. 

Electro-mechanics 

Our model for cardiac contractile mechanics is based on the approximate model of 
cooperative activation and crossbridge cycling reported by Rice et al. [55] with the fol- 
lowing modifications: (a) the first-order rate constants for the transformation of the 
troponin/tropomyosin regulatory complex (outside the single overlap region between 
the thick and thin filaments) from a crossbridge non-permitting state to a crossbridge 
permitting state and vice-versa are chosen as 500 s-1 and 50 s-1 respectively in order 
to reproduce results reported by Rice et al. [55]; (b) the ^-adrenergic agonist isopro- 
terenol (ISO) is known to cause a decrease in myofilament Ca^^ sensitivity as a result 
of PKA mediated phosphorylation of troponin I [51,56] at Ser23/Ser24. Specifically, a 
two-state Markovian model is added to allow cAMP-dependent PKA-mediated inter- 
action between troponin I (Tnl) and the 0/3^"^+ -binding regulatory site on troponin. As 
reported by Messer et al. [57], the unphosphorylated form of Tnl (Tnl^^) modulates 
the Ca?'^ affinity of the regulatory site on troponin. We model the effects of cAMP 
by allowing the cumulative activation rate constant for Ca^^-hmdimg to the troponin 
regulatory site to be a function of unphosphorylated Tnl^^, the availability of which is 
in turn dependent on the amount of [cAMP] present (Appendix, Equations 10,12); (c) 
the large Qio values used by Rice et al. {Qiappy Qh/, Qh^ and Qg;,^^, Table 1, [55]) are 
decreased from 6.25 to 2.25 in order to reproduce temperature dependence of peak force 
developed in intact thin rat ventricular trabeculae [58]. The rate constants (Appendix, 
Equations 13,17) governing cross-bridge kinetics are modeled as functions of [cAMP] to 
reproduce stimulation frequency dependent increase in contraction and relaxation rates. 
Although a calmodulin (CaM) mediated pathway has been reported [59] to be responsi- 
ble for modulation of myofibrillar C<3^^+ -sensitivity (implying a possible CaM mediated 
role for Ca-dependent kinases or phosphatases in regulating myofilament contractility, 
particularly in frequency dependent acceleration of relaxation), we refrain from modeling 
this effect as the molecular mechanisms involved remain unresolved. 

Results 

From our modeling standpoint, the dyadic coupling unit (DCU) as defined by Krishna 
et al. [15] is a fundamental element involved in the mechanism of CICR. These previous 
studies have described the control features of this unit, as well as its interaction with the 
SERCA pump and free sarcolemmal pumps and exchangers to achieve a homeostatic reg- 
ulation of myoplasmic Ca^^ concentration. We now extend our voltage clamp studies to 
address the subject of frequency-dependent characteristics of CICR and begin with the 
study of frequency dependence of the DCU, one of the most important components of the 
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model All the frequency-dependent behavior discussed here is at steady state (150 cycles) 
unless otherwise specified. The first task is to examine the frequency dependence of the 
^Ca,L trigger current, followed by the RyR Ca^^ channel, highlighting the rate-dependent 
CaM-mediated signalling involved. This analysis is followed by examining the SR C^2+ 
content and the various factors controlling it in a rate-dependent manner. A quantitative 
study of the overall cellular Ca^^ balance is performed to highlight its rate-dependent 
feature. Emphasis is placed on relative roles played by the longitudinal sarcoplasmic retic- 
ulum (LSR) membrane SERCA pump and the plasma membrane Na^ ICa^^ exchanger, 
as two principal Ca^^ transport routes in the maintenance of Ca^^ homeostasis. We 
finally examine the myoplasmic Ca^^ transient as a function of frequency with a particu- 
lar interest in the rate dependence of the force-frequency response (FFR) generated by the 
coupled electromechanical model. This is subsequently followed by an investigation of 
the rate-dependent influence of cAMP-mediated -adrenergic stimulation on the cardiac 
contractile response. 

L-type Ca^+ current (lca,L) 

An increase in stimulation frequency from 0.5 Hz to 8 Hz, results in a frequency- 
dependent monotonic increase in the peak trigger current (inset in Figure 4A and 
Figure 4B) while slowing down the rate of decline after it reaches its maximum (note 
the traces corresponding to 0.5 Hz to 8 Hz in Figure 4A). The most critical mechanism 
involved in frequency encoding of the Ica,L channel activity is the rate-dependent change 
in the average level of activated-CaMKII (Figure 4C), which is known to assist CaM- 
mediated C<3^^+ -dependent facilitation (CDF). As stimulation frequency is increased from 
0.5 Hz to 8 Hz, the increase in peak /c^,/^ (Figure 4B) closely tracks the increase in the 




Time (ms) Frequency (Hz) 



Figure 4 Ica.L- Frequency-dependent characteristics of tine DHP-sensitive Co^^ cliannel Ica.L, wliicli forms tine 
trigger current. Voltage clamp protocol used is a 50 ms step pulse to 10 mVfrom a holding potential of -40 
mV. (A) Icai at different frequencies of stimulation ranging from 0.5 Hz to 8.0 Hz. Inset gives an expanded 
view of the rate-dependent modulation oUcai over the range of frequencies investigated (B) Dependence 
of peak lca,L current on frequency of stimulation (C) Frequency-dependent increase in average levels of 
activated auxiliary proteins CaMKII and CaN in the dyad (D) Frequency-dependent modulation of peak open 
probability of the DHP-sensitive Cc/^+ channel. (E) Rate-dependent changes in Icai upstroke velocity. 
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average level of activated-CaMKIL However, beyond 8 Hz a decrease in peak is observed 
(Figure 4B) despite a further increase in CaMKII (Figure 4C). This occurs as a result of 
incomplete channel recovery at high (> 8Hz) stimulation rates, which results in a decline 
in peak channel open probability (Figure 4D). At low (0.5 Hz < f < 4.0 Hz) stimulation 
rates, an increase in activated CaN (85% to 97%) is also known to enhance Ica,L chan- 
nel activity, whereas at higher (> 4 Hz) rates, the lack of a substantial rate-dependent 
increase in its average level (Figure 4C) minimizes its role in C<3J^+ -dependent facilitation. 
The maximum value attained by the open probability of the DHP-sensitive Ca^^ chan- 
nel (Figure 4D) reflects the trend shown by the peak value of the trigger current over the 
entire range of stimulation frequencies (0.5 to 12.0 Hz) investigated. At frequencies less 
than 8 Hz, the upstroke velocity oiIca,L current can be seen to increase with an increase in 
frequency, but above 8 Hz it begins to decline (Figure 4E) due to insufficient time for full 
channel recovery. It is important to note that the model predicts a frequency-dependent 
modulation in peak Ica,L current of less than 20% over the entire frequency range (0.5 Hz 
to 12.0 Hz). This small modulation of peak current (Figure 4B) is far less than the per- 
centile changes in CaMKII activation (50%, Figure 4C), due to the insufficient time for 
channel recovery at high (> 4 Hz) stimulation rates. 

RyR Ca^"*" release 

Figure 3A shows a phase plot of RyR open probability {Pq) versus [Ca^^]KyR constructed 
from model-generated data corresponding to different stimulation rates with the inset 
showing the peak open probability attained by each of the phase loops. As the stimulation 
frequency is increased from 0.5 Hz to 4 Hz, a marginal increase (< 1%) in peak RyR open 
probability occurs (inset in Figure 3A) as a result of the factors: (a) a frequency-dependent 
CaMKII mediated (Figure 3B) Ica,L facilitation (Figure 4A); (b) a moderate increase in SR 
Ca?'^ content (despite a relatively constant maximal SERCA uptake rate), due to increased 
trigger current and the subsequent increase in [Ca^^]myo combined with an increase in 
[Na^]myo (30% in Figure 3B) which impedes Ca^^ extrusion via the NCX; and (c) a direct 
CaMKII-mediated enhancement in RyR release. As the stimulation frequency is increased 
from 4 Hz to 8 Hz, despite a moderate increase in both the trigger current and SR Ca^^ 
content, a small decrease in RyR open probability occurs as a result of an increase in 
post-release [Ca?''^]jSR (Figure 3C), which forces an incomplete luminal sensor-based RyR 
recovery (as described in Krishna et al. [15]). Beyond 8.0 Hz, the peak RyR open probabil- 
ity decreases due to two mechanisms: (a) a small parallel decrease in the trigger current 
(Figure 4B) indicating a close coupling enforced by a stable CICR, and (b) insufficient 
time for full channel recovery, accompanied by a falling pre-release diastolic jSR CcP'^ 
level (Figure 3C). Both the declining trigger current and the declining [Ca^^]jSR result in 
a strong decline in [Ca^^]RyR, The decrease in [Ca^^]RyR at high (> 8 Hz) frequencies can 
be seen in Figure 3A, where the open probability loop for 10 Hz is enclosed within that of 
8 Hz. The frequency dependence of the area enclosed by the loops (which indicates the 
amount of SR Ca^^ released into the dyad) mirrors that of the peak RyR open probability. 
The SR Ca^^ release from the Ry-sensitive receptor exhibits strong frequency-dependent 
behavior as shown in Figure 3D. Instantaneous RyR flux is obtained by multiplying the 
open probability by the concentration gradient across the channel. The rate dependence 
of the peak instantaneous RyR flux (Figure 3D, E) is a result of two factors: (a) the fre- 
quency dependence of peak RyR open probability (Fig 3A), which indicates the degree of 
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recruitment of RyR release channels; and (b) the pre-release SR Ca?'^ content (Figure 3C), 
which establishes the initial concentration gradient across the release channel. Figure 3D 
also shows the gradation in the time required for RyR recovery. 

The frequency dependence of RyR release activation (Figure 3F), shares its character- 
istics with the peak RyR open probability (Figure 3A). Besides increasing trigger current, 
RyR release is also facilitated by increasing levels of activate d-CaMKII (Figure 3E) and 
increasing SR Ca^^ content (Figure 3C). The frequency-dependent modulation of RyR 
inactivation rate (Figure 3F) mimics the rate dependence of minimum SR Ca}^ levels 
reached after release, mediated via the luminal sensor. With increasing stimulation fre- 
quency (0.5 Hz < f < 8.0 Hz), despite decreasing time available for uptake, increasing 
pre-release SR Ca^^ content is achieved only by increased rate of SR filling. This increased 
rate of SR filling translates into a faster decline in inhibition via the luminal sensor 
(a RyR inactivation mechanism described in Krishna et al. [15]). Hence, the rate of inacti- 
vation decreases as the frequency is increased from 0.5 Hz to 8.0 Hz (Figure 3F) owing to 
increasing SR Ca^^ content in this range of frequencies. With increase in frequency, the 
rate-dependent increase in the level of activated-CaMKII in the dyadic space also causes 
an increased CaMKII-mediated upregulation of the ryanodine receptor [40,41,43,45-48]. 
This CaMKII mediated up-regulation delays the onset of declining SR CcP'^ content- 
driven (8.0 Hz < f < 12.0 Hz) increase in the rate of RyR channel inactivation. The 
RyR channel experiences a frequency-dependent modulation by both the trigger current 
and the amount of activated-CaMKII on its dyadic side. On its luminal side, it experi- 
ences modulation by the luminal sensor (as described in Krishna et al. [15]) controlling 
refractoriness, and the SR Ca^^ content providing the drive for Ca}'^ through the channel. 

SR Ca2+ content 

The frequency dependence of pre-release Ca?'^ level in the SR is a result of several 
factors namely: (a) available [Ca^'^]myo for sequestration; (b) the rate-dependent behav- 
ior of the SERCA pump; (c) the frequency dependence of the release characteristics of 
Ry-sensitive Ca?'^ channel; and (d) the cumulative transmembrane Ca^^ transport via 
ica,L and NCX. The combined effect of these factors results in the biphasic relationship 
(inset in Figure 5A) between free Ca^^ content in the SR and the frequency of stimula- 
tion. As the frequency is increased from 0.5 Hz to 8.0 Hz, the SR Ca^^ content increases 
due to the following frequency-dependent, active CaMKII-mediated effects: (a) enhance- 
ment of maximal uptake rate of the SERCA pump due to increase in PLB phosphorylation 
(Figure 5B) assisted by increasing levels of activated-CaMKII (Figure 5C); (b) decrease in 
the half-activation constant for the forward (myoplasm to SR) operation of the SERCA 
pump (Figure 5D) translating into increased Ca^^ sensitivity of the pump for uptake; and 
(c) increase in half-activation constant for the backward (SR to myoplasm) operation of 
the SERCA pump (Figure 5D), reducing tendency for back-flow via the SERCA pump. 
The increase in SR Ca^^ content is also facilitated by a CaMKII mediated enhancement 
in Ica,L as well as inhibition of Ca?'^ extrusion via NCX due to increasing [Na^]ynyo' As 
the stimulation frequency is gradually increased further from 8.0 Hz to 12.0 Hz, a steep 
decrease in SR Ca^^ content is observed (inset in Figure 5A). This frequency-dependent 
decrease in pre-release SR Ca^^ content at very high (> 8.0 Hz) stimulation frequen- 
cies is a result of the following: (a) a decrease in maximal uptake rate of the SERCA 
pump (Figure 5B) along with a significant decrease in time available for resequestration 
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Figure 5 [Ca^+lysff.Mechanisms underlying the frequency-dependent modulation in SR Cci^+ content. 
Voltage clamp protocol used is a 50 ms step pulse to 10 mVfrom a holding potential of -40 mV. (A) [Ca^'^]jSR 
traces indicating rate-dependent changes in SR filling. The inset shows the frequency dependence of the 
maximum and minimum SR Cc/^+ concentration. (B) Frequency dependence of the maximal pumping rate 
of SERCA and the average level of phosphorylated PLB. (C) Frequency-dependent increase in average levels 
of activated auxiliary proteins CaMKII and CaN in the myoplasm. (D) Frequency dependence of the 
half-activation constant for the forward and backward flow via the SERCA pump. (E) Rate dependence of the 
relative role of the SERCA pump (when compared to the No'^/Ca^^ exchanger (//voCo)and the plasma 
membrane Ca^^ ATPase pump {Ipmca)- 



of cytosolic Ca^^; (b) a decrease in Ca^^ entry into the cell via the trigger current Ica,L 
(Figure 4B); and (c) relatively constant intracellular [Na^]myo (Figure 3B). Between stim- 
ulation frequencies 0.5 Hz and 8.0 Hz, the characteristics of the post-release SR Ca^^ 
level tracks the pre-release peak SR Ca^^ content. However, beyond 8.0 Hz, a declining 
SR Ca^^ content causes increasing inhibition on RyR release via the luminal sensor [15], 
resulting in a gradual increase in post-release SR Ca^^ level. 

The frequency encoding involved in the uptake mechanism via the SERCA pump is 
strongly regulated by two key proteins, PLB and activated-CaMKII. The rate-dependent 
decrease in the level of unphosphorylated phospholamban (PLB) (Figure 5B) capable 
of inhibiting the SERCA pump is a result of increasing levels of activated-CaMKII 
(Figure 5C). It is important to note that phosphorylation of PLB is sensitive to small 
changes in activated-CaMKII [23]. Our results show that activated-CaMKII (Figure 5C) 
has a prominent role in rate-dependent modulation of the SERCA pump despite its low 
level (three orders of magnitude compared to activated CaN (note the 10~^ on left ordi- 
nate in Figure 5C)) in the myoplasm [2]. Our results show that activated CaN (Figure 5C) 
has minimal role in rate-dependent inhibition of the SERCA pump. However, over- 
expression of CaN in failing heart is known to significantly compromise SERCA activity 
[54]. At a given frequency of stimulation, an increase in trigger current forces a grad- 
ual increase in SR Ca^^ uptake which in turn results in a corresponding increase in 
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RyR release, which translates into an elevated cytosolic Ca^^ level and subsequently an 
increased activated-CaMKII level. This causes an increase in SR Ca^^ uptake rate via 
phosphorylation of PLB, thus completing a positive feedback loop. Hence, the peak Ica,L 
trigger current sets the peak SR Ca^'^ content if the SERCA pump is not operating in satu- 
ration. Thus the frequency dependence of the maximal rate of uptake by the SERCA pump 
(Eigure 5B) shares its characteristics with peak Ica,L current (Eigure 4B) linked together 
by activated-CaMKII (Eigure 5C). It is important to note that rate-dependent enhance- 
ment in SERCA activty by CaMKII significantly overrides the inhibitory influence of CaN 
in a non-failing heart [54]. Eigure 5B also shows the underlying rate-dependent decline in 
unphosphorylated PLB responsible for rate-dependent increase in SR uptake. 

Although marginally (< 1% change over 0.5 Hz < f < 12.0 Hz) frequency-dependent, the 
SERCA pump has a primary role in resequestration of myoplasmic Ca^^ when compared 
with the combined influence oiNa^ ICa?'^ exchange (InuCo) as well as Ca^~^ transport via 
the plasma membrane Ca^^ ATPase pump, (/pmca)- As the frequency of stimulation is 
increased from 0.5 Hz to 4.0 Hz, increasing activated-CaMKII-dependent enhancement 
in SERCA activity assisted by a decrease in Ca^~^ extrusion via NCX due to increasing 
[Na'^]ynyo (Eigure 3B) result in a marginal increase in the relative role of SERCA pump 
as shown in Eigure 5E. However, as frequency is further increased, declining SERCA- 
mediated uptake activity due to a decrease in CaMKII-mediated SERCA enhancement, as 
well as reduced time available for uptake results in a small negative slope to the frequency 
dependence of the relative role of SERCA pump (Eigure 5E). 

Myoplasmic Ca^+ transient [Ca^+]myo 

Eigure 6 A shows the traces of [Ca^'^]niyo at increasing frequencies showing an increase 
in the peak myoplasmic Ca^^ transient during systole with increasing stimulation fre- 
quency in the range 0.5 Hz to 8.0 Hz. A 'primary phase' negative peak EER (a decrease in 
peak contractile force with increase in stimulation frequency observed at low (< 1.0 Hz) 
stimulation rates) reported in studies on rat ventricular myocytes [8,9,60,61] could be 
a result of low unphysiological temperatures employed in those studies. The inset in 
Eigure 6A shows that the peak rate of rise {Rrise) and decay i^decay) of the [Ca?'^]myo tran- 
sient increases with the frequency of stimulation from 0.5 Hz to 8.0 Hz. The increase in 
rate of rise is a result of an increase in peak trigger current Ica,L (Eigure 4B) and SR Ca^^ 
content (inset in Eigure 5A) while the increase in rate of decay is a result of CaMKII- 
mediated acceleration of relaxation (Eigure 5C) due to a rate-dependent enhancement in 
uptake by the SERCA pump. An opposite effect is observed in the rate of rise and decay 
of the [Ca^'^]myo transient beyond a stimulation rate of 8.0 Hz due to declining peak trig- 
ger current and reduced CaMKII activity. As shown in Eigure 6B, increasing stimulation 
rate from 0.5 Hz to 8.0 Hz results in a steep increase in peak [Ca^~^]myo (maximum), 
which parallels a strong increase in pre-release diastolic [Ca^^]jSR (Eigure 5A). This 
increase in [Ca^^]jsR is mirrored by a corresponding increase in pre-release minimum 
value attained by [Ca^'^]myo during diastole (Eigure 6B). As the stimulation frequency is 
increased beyond 8.0 Hz, although the SERCA activity is saturated [12], the declining SR 
Ca^^ content (inset in Eigure 5 A) aids in the decrease of peak [Ca^~^]myo by compromising 
SR release. However, a lack of sufficient time for uptake continues to manifest in increas- 
ing pre-release diastolic \C(P'^\yyiyo (Eigure 6B) at these rates. Our study suggests that a 
decrease in pH/ causing a decrease in myofilament Ca}^ sensitivity [62] is not required 
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Figure 6 [Ca^+lmyo- Frequency-dependent modulation of myoplasmic C6/^+ transient. Voltage clamp 
protocol used is a 50 ms step pulse to 10 mVfrom a holding potential of -40 mV. (A) [Cci^+] inyQ traces 
showing rate-dependent effects on systolic Ca^+ concentration in the myoplasm for increasing stimulation 
frequencies from 0.5 Hz to 8.0 Hz. The inset shows the change in the maximum rate of rise and decay of the 
[Ca^+]myo transient with stimulation frequency. (B) Frequency dependence of the maximum and minimum 
myoplasmic Ca^+ concentration. (C) Bi-phasic rate dependence of EC coupling gain. 



to cause the secondary phase' negative peak FFR (a decrease in peak contractile force 
with increase in stimulation frequency observed at high (> 8.0 Hz) stimulation rates). 
We attribute this phase which involves a decrease in force response, to the declining 
peak \Ca}^\ynyo as seen in Figure 6B. The rate-dependence of peak \Ca}^\ynyo for increas- 
ing stimulation frequencies (0.5 Hz < f < 12.0 Hz) is comparable to that reported by 
Kentish et al. (Figure 1, [63]) in rat ventricular trabeculae. Figure 6C shows the frequency- 
dependent changes in gain calculated as the ratio of peak Ca}^ transient in the presence 
of CICR to the peak calcium transient in its absence, contributed by the trigger calcium 
alone [64]. The EC coupling gain increases with the frequency of stimulation from 0.5 
Hz to 8.0 Hz, due to rate dependent CaMKII-mediated increase in uptake rate of the 
SERCA pump (Figure 5B), which decreases the peak [Ca^~^]myo in the absence of CICR 
and increases peak [Ca^~^]myo in its presence (due to increase in SR content (inset in 
Figure 5A)). Beyond a stimulation frequency of 8.0 Hz, the rapid decrease in SR content 
causes a steep decline in EC coupling gain. 

Positive force-frequency relationship 

Changing the stimulation rate has a very prominent effect on cardiac muscle contrac- 
tion by way of affecting the amplitude and time course of the intracellular Ca^^ transient 
(Figures 6 and 7A). Here, we simulate a positive peak force-frequency response elicited 
using a voltage clamp stimulation protocol. Figure 7B shows the traces for the isometric 
contractile force generated for increasing stimulation frequencies. As stimulation fre- 
quency increases from 0.5 Hz to 8.0 Hz, peak contractile force increases (Figures 7B 
and C) mirroring an increase in peak [Ca^~^]myo' As seen in Figures 4B and 5B, this is 
mainly a result of CaMKII-dependent upregulation of Ica,L and SERCA activity which 
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Figure 7 Contractile force and sarcomere length. Rate-dependent changes in myoplasmic Cci^+ 
concentration, tine corresponding contractile force generated and tine resulting sarcomere length. Voltage 
clamp protocol used is a 50 ms step pulse to 10 mVfrom a holding potential of -40 mV. (A) Frequency 
dependence of the maximum and minimum myoplasmic Cc/^+ concentration (B) Traces for normalized 
isometric force developed at stimulations rates of 4, 5, 6, and 8 Hz. (C) Rate dependence of the maximum and 
minimum normalized force generated. (D) Phase loops for normalized isometric force vs [Ca^+]myo for 
increasing frequency of stimulation (4.0 Hz to 8.0 Hz) corresponding to panel B. (E) Traces for sarcomere 
contraction corresponding to isometric force in panel B. (F) Frequency dependence of minimum sarcomere 
length achieved at peak contraction. The legend for traces in panel B is shared by panel D and E. 



translates into a frequency-dependent increase in SR Ca^^ content, and thereby SR 
Ca^^ release. As reported in the literature [55], the non-linear sigmoidal steady state 
Force- Ca^^ relationship (overlayed in Figure 7D) causes a non-linear response in peak 
force developed for a corresponding change in peak [Ca^~^]myo' As shown in Figure 7D the 
contraction-relaxation coupling point (CRCP) moves to larger values of force and Ca^^ 
concentration with increase in rate of stimulation. The peak contractile force is attained 
at a physiological rate of 8.0 Hz, beyond which the gradual decline in peak Ica,L cur- 
rent (Figure 4B), which forces a corresponding decrease in SR Ca^^ content (Figure 5A), 
and hence RyR release (Figure 3E), causes a marginal decline in peak contractile force 
generated. The frequency at which peak contractile force is generated is temperature- 
dependent as reported by Layland et al. [12] and shifts to the right on the frequency axis 
as the temperature increases. The peak force frequency response obtained mimics the 
experimentally observed frequency dependence (Figure 1, [63]). The minimum value of 
contractile force per cycle (end-diastolic) increases with increasing rate of stimulation 
as a result of the combination of factors: (a) increasing myoplasmic Ca^^ levels due to 
rate-dependent increase in SR release (Fig 6); and (b) decreasing time available for rese- 
questering released Ca^^ from the myoplasm. The frequency dependence of the force 
response translates into a reciprocal rate dependent characteristic of the sarcomere length 
as shown in Figures 7E and F. 
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Effect of maximal jS-adrenergic stimulation 

y6 -adrenergic stimulation bears a very crucial physiological relevance in normal function- 
ing of the heart. We have investigated the effect of ^-adrenergic stimulation on [Ca^^]ynyo 
under voltage clamp conditions. The degree of CaMKII mediated Ica,L facilitation is 
kept identical to the control case discussed earlier. Besides a moderate CaMKII-mediated 
[23,65] enhancement in Ica,u we observe that increase in stimulation frequency causes a 
strong -adrenergic stimulation-dependent facilitation in Ica,Ly resulting in an enhanced 
positive slope in the FFR as the stimulation frequency is increased. These results tend 
to agree with studies that show that increased cAMP levels reverse the abnormalities 
in the FFR that are found in end-stage heart failure [66], and that moderate stimula- 
tion by isoproterenol partly reverses the negative FFR found in NYHA (New York Heart 
Association) class IV (severe) type heart failure and preserves the positive FFR in non- 
failing myocardium [67]. It is important to note that variations in experimental protocol 
could result in undesirable artifacts owing to the time-sensitive nature of the effects of 
cAMP [68]. Here, we show that the presence of maximal ^-adrenergic stimulation results 
in a stronger positive FFR (frequency dependence of peak contractile force generated) at 
physiological stimulation frequencies. 

At low (< 2.5 Hz) stimulation frequencies ^-adrenergic stimulation plays no role in 
altering intracellular Ca^^ dynamics. However, in the presence of maximal ^-adrenergic 
stimulation (elevated levels of cAMP), an increase in stimulation rate (> 2.5 Hz) causes 
increased up-regulation oiIca,L [69-71]. The model predicts (Figure 8A) that as the fre- 
quency of stimulation is increased in steps, peak Ica,L magnitude increases at a regular 
rate, before declining at high stimulation frequencies (> 8.0 Hz) as a result of incom- 
plete channel recovery. As shown in Figure 8B, as the frequency is increased from 2.5 
Hz to 8.0 Hz, the peak SR Ca^^ content increases due to two mechanisms, namely 
cAMP-dependent enhancement of both trigger current and a rate-dependent increase 
in uptake by the SERCA pump as a result of a stronger decline in the average level of 
unphosphorylated PLB, which stimulates the uptake. It is important to note that although 
the increase in peak [Ca^~^]myo (traces in Figure 8C) causes a strong CaM-mediated 
CaMKII-dependent up-regulation of the SERCA pump due to an increase in PLB phos- 
phorylation, cAMP-mediated phosphorylation of PLB causes further enhancement of the 
uptake mechanism. An increase in stimulation rate beyond 8.0 Hz, causes a decrease in 
SR Ca^^ content due to insufficient time for SR recovery along with a saturated operation 
of the SERCA pump. The presence of ^-adrenergic stimulation significantly alters the 
relationship between the peak force response and the frequency of stimulation. As shown 
in the traces for force response in Figure 9A, besides an increase in peak force gener- 
ated (a result of increased peak \Ca^^\myo (Figures 8C and 9B)), ^-adrenergic stimulation 
causes an increase in contraction and relaxation rate. This cAMP-mediated increase in 
the rate of sarcomere contraction and relaxation is a combined result of increasing rate of 
rise and decline of the C^s^"^"^ -transient and a cAMP-mediated enhancement in rate kinet- 
ics governing cross-bridge formation. The inset in Figure 9A is a plot of the maximum 
rate of relaxation i^relaoc) versus maximum rate of contraction ^act) for increasing stim- 
ulation frequency with and without ^-adrenergic stimulation (shows significant overlap). 
This linear relationship highlights contraction-relaxation coupling, and represents a key 
intrinsic property of the contractile myofilaments [72]. As shown in Figure 9B, although 
^-adrenergic stimulation further assists (due to the upregulation of lca,L and SERCA 
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Figure 8 Effect of maximal j8-adrenergic stimulation on Ca^'^ dynamics. A comparison between 
cAMP-mediated maximal yS-adrenergic stimulation and basal condition. Voltage clamp protocol used is a 50 
ms step pulse to 10 mVfrom a holding potential of -40 mV. (A) Rate dependence of peak/co,/_ current. (B) 
Frequency dependence of the maximum pre-release SR C6/^+ concentration. (C) [Ca^+]myo traces showing 
rate-dependent effects on systolic Ca^+ concentration in the myoplasm for increasing stimulation 
frequencies from 4.0 Hz to 8.0 Hz. 



pump) a CaMKII-mediated rate -dependent increase in peak [Ca'^~^]^yo, CaMKII is the 
key factor responsible for this behavior (a 2-fold increase in peak [Ca^^]myo from 4 Hz 
to 8 Hz both in the presence/absence of cAMP-mediated effects). In the presence of 
-adrenergic stimulation, Figure 9C shows an enhancement in peak force generated and 
a decrease in minimum force as a result of a parallel trend in the peak systolic and end- 
diastolic levels of [Ca^~^]myo (Figure 9B) supported by faster cross-bridge kinetics. The 
delta increment in force due to ^-adrenergic stimulation (compared to the control case) is 
maximum at the stimulation rate of 5 Hz (Figure 9C) and decreases both above and below 
this frequency. The cAMP-mediated effect on the dependence of maximum sarcomere 
contraction (Figure 9D) on stimulation frequency mirrors that of the peak force response 
(Figure 9C) showing an opposite trend. 

Discussion 

In this study, we have examined the role of two key mechanisms: (a) CaMKII-mediated 
upregulation; and (b) frequency-dependent cAMP-mediated ^-adrenergic stimulation 
on Ca^^ cycling and their effect on the force-frequency response. In particular, we 
have modeled their two crucial signalling mechanisms: the rate-dependent upregula- 
tion of the DHP-sensitive Ica,L channel [25,69,73,74] and SERCA pump activity [75,76], 
which our study suggests are the key factors responsible for the characteristics of FFR 
in rat ventricular myocytes. These two signalling pathways mediate Ica,L enhancement 
bringing an extra supply of trigger current into the dyad to both enhance CICR and 
assist in a rate-dependent build up in pre-release diastolic SR Ca^^ concentration. This 
frequency-dependent increase in SR Ca^^ content is strongly assisted by a -adrenergic 
stimulation-dependent cAMP-mediated increase in PLB phosphorylation, which relieves 
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Figure 9 Effect of maximal jS-adrenergic stimulation on force response. A comparison between 
cAMP-mediated maximal yS-adrenergic stimulation and basal condition. Voltage clamp protocol used is a 50 
ms step pulse to 10 mVfrom a holding potential of -40 mV. (A) Traces for normalized isometric force 
developed for increasing stimulation frequencies from 4.0 Hz to 8.0 Hz. The inset shows a plot of the 
maximum rate of relaxation (Rreiax) versus maximum rate of contraction (Ract) fot" increasing stimulation 
frequency. (B) The frequency dependence of the maximum and minimum myoplasmic Cc/^+ concentration. 
(C) Rate dependence of the maximum and minimum normalized force generated. (D) Frequency 
dependence of minimum sarcomere length achieved at peak contraction. 



the inhibition exerted on the SERCA pump. This is especially important in the light of 
decreasing time for Ca^^ uptake with increasing stimulation rate. 

We have systematically analyzed the important role of C<3^^+ -dependent CaM-mediated 
function of CaMKII and CaN as rate sensors. CaMKII is known to regulate EC coupling 
[41,45] and has substrates in both the dyadic domain {Ica,L and Iryr) [32,41] as well as 
the myoplasm (phospholamban) [77], and immunoprecipitates with the Ica,L channel [32] 
and RyR [78]. Our study confirms that activated-CaMKII causes rate-dependent accel- 
eration of relaxation [3] both by direct phosphorylation of the SERCA pump as well as 
phosphorylation of PLB resulting in a decreased inhibition of the SERCA pump. Although 
the time course of the indirect effect of CaMKII via PLB phosphorylation is reported 
to be delayed with regard to frequency-dependent acceleration of relaxation [79], the 
direct phosphorylation of the SERCA pump by activated-CaMKII could cause an imme- 
diate increase in maximal pumping rate [4,23]. Despite poor sensitivity of CaMKII to 
myoplasmic Ca^^ signals due to its relatively low affinity for CaM and an inadequate sup- 
ply of myoplasmic Co^^-CbM complex (compared to the dyad), activated-CaMKII in the 
cytosol plays an important role in the enhancement of uptake via the SERCA pump except 
at high (> 8 Hz, inset in Figure 5A) stimulation frequencies. Here, the rapid decrease 
in time available for uptake results in gradually declining pre-release SR Ca^^ content. 
CaMKII has a more critical role in the dyadic cleft where we show [2] substantial CaMKII 
activity owing to large local Ca^^ concentration and CaM enrichment. It is reported that 
frequency-dependent CaMKII activity contributes to Ica,L channel facilitation [31] and 
assists a positive force-frequency relationship [80]. Our study shows that CaMKII causes 
a significant enhancement oiIca,L even under basal conditions (absence of ^-adrenergic 
stimulation) as the stimulation frequency is increased to 8.0 Hz (Figure 4B), beyond which 
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incomplete recovery begins to cause a decline in peak Ica,L' A rate-dependent increase in 
SR Ca^^ content can be achieved only when the CaMKII-mediated (Figure 4C) increase 
in trigger current (Figure 4B) is supported by an increase in uptake via the SERCA pump 
by activated-CaMKII in the myoplasm (Figure 5C) which enhances PLB phosphorylation. 
CaMKII fails to sustain a positive FFR at very high (> 8.0 Hz) stimulation frequencies due 
to incomplete channel recovery causing a decline in trigger current and subsequently SR 
Ca^^ content. Data from atrial myocytes [81] show that myoplasmic CaN is also found to 
be responsive to increased pacing frequency. 

Although CaN is also involved in Ica,L enhancement in the dyad at low (< 2.0 Hz) stim- 
ulation frequencies, at high stimulation rates (> 2.0 Hz), the large Ca^^ signals along with 
very slow dissociation of the Ca^CaM-C^N complex force it to be constitutively active, 
limiting its modulatory role. The presence of y6 -adrenergic stimulation, which also sup- 
ports a sustained frequency-dependent increase in Ica,L (Figure 8A) and SERCA activity 
(Figure 8B) results in enhanced contractile response (Figure 9 A and C). 

Here, we have integrated an electro-mechanical model with our electrophysiological 
model of a rat ventricular myocyte [15] to describe isometric contractile force gener- 
ated by myofibrils and investigate its rate-dependent characteristics. Our model for the 
cardiac contractile mechanism is derived from the description of cooperative activation 
and cross-bridge cycling given by Rice et al. [55]. However, we have incorporated a rate- 
dependent cAMP mediated enhancement in cross-bridge kinetics which is shown to be 
responsible for sustaining linear contraction-relaxation coupling (inset in Figure 9A) at 
increasing frequencies. It is important to note that the presence of CaMKII-mediated 
^Ca,L enhancement or ^-adrenergic stimulation results in physiologically relevant rate- 
dependent increase in peak cardiac contractile force (positive FFR). 

Model limitations 

1. This model of a rat ventricular myocyte is limited to Ca^^ related channel, 
exchanger and pumps {Ica,L^ ^NaCa^ hMCAi Uyr and SERCA pump), while lacking 
exclusive Na'^ or related channels and transporters. Hence, it is aimed at 
mimicking voltage clamp conditions where channels other than calcium are 
blocked, and it cannot be used to study any action potential-induced Ca^~^ 
transient-frequency relationship. However, its focus on the Ca^'^ dynamics alone 
allows one to comprehend more clearly the important role of various 

C<3^^+ -dependent regulatory proteins such as CaMKII, CaN, PLB and cAMP in 
affecting multiple targets and thus generating a cell's response to a change in the 
frequency of stimulation. 

2. It is known that CaMKII alters the function of numerous ion channels and Ca^^ 
regulatory targets in a rate-dependent fashion. However, disparate findings exist 
on its modulation of targets such as the SERCA [82,83] and RyR channel 
[40,44,45]. Although our model aligns with rate-dependent CaN-mediated 
inhibition of the SERCA pump, its role in modulating SERCA activity is 
controversial due to conflicting findings [53,54]. Similarly, the effect of 

-adrenergic stimulation via cAMP particularly, its dose-dependent influence on 
L-Type Ca^^ channels both in terms of modifying the single channel behavior 
such as Ca^'^ ion permeability (channel conductance) as well as overall channel 
recruitment characteristics (open probability) is not clearly understood. While 
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cAMP has been shown to increase channel open probability [35,36], increased 
levels of cAMP also result in increased phosphorylation of L-type Ca^^ channels, 
causing an increased permeability to Ca^^ ions [36-38]. Although we have 
modeled enhancement in Ica,L due to beta-adrenergic stimulation solely by a 
rate-dependent cAMP-mediated increase in channel conductance, the relative 
contribution of these factors to cAMP-dependent Ica,L enhancement has to be 
examined further. A detailed investigation is required to clarify the nature of these 
interactions. We hope that this study would help motivate more pointed 
experimental investigation of frequency-dependent CaMKII, CaN and cAMP 
effects on FFR in rat ventricular myocytes. 

3. The effect of -adrenergic stimulation on cardiac Na'^ ICa^^ exchange has been 
controversial. Perchenet et al. [84] report an enhancement in Na^ I Ca^^ exchange 
by the yS-adrenergic/PKA-mediated phosphorylation of the exchanger protein. 
However, a cAMP-mediated enhancement in NCX activity would impede rate 
dependent increase in SR Ca^^ content. We base our model for the Na^ ICa^^ 
exchanger on a more recent study by Lin et al. [85] which supports the view that 
^-adrenergic stimulation does not upregulate Na^ ICcP'^ exchange current. 

4. The cooperative activation of the thin filament and the strain-dependent 
transitions of the cross-bridge cycle have been approximately modeled as 
non-spatial, state-variables. However, this simplification is valid as these 
transitions are inherently local phenomena and the model reproduces a wide 
range of steady state and dynamic responses in cardiac muscle [55]. 

Conclusion 

Using a mathematical model of an isolated rat ventricular myocyte in a voltage clamp 
setting, we have systematically examined the issue of rate-dependence in the proper func- 
tioning of the dyadic coupling unit, the regulation of SERC A function to provide adequate 
SR Ca^^ content, the peak amplitude of the myoplasmic Ca^^ transient and the com- 
plex interaction of all these factors. Given the complexity of these interacting systems, 
computer modeling gives an insight into the relative roles of different Ca^^ transport 
mechanisms. Our simulations explain the 0/3^"^+ -dependent, CaM-mediated, rate sensitive 
effects of CaMKII and CaN on various intracellular targets. We also investigate a signif- 
icant, frequency-dependent, cAMP-mediated effect of -adrenergic stimulation and its 
modulatory influence on the Ica,L channel as well as the SERCA pump. Rate-dependent 
CaMKII mediated Ica,L facilitation as well as cAMP-dependent upregulation of intracel- 
lular targets could play a vital role in reversing the negative FER found in failing hearts. 
However, further studies are required to develop a clear understanding of the relative 
role of CaMKII and cAMP in the rate-dependent up-regulation of various intra-cellular 
targets especially the DHP-sensitive Ica,L channel. This would help in assigning rate- 
dependent weights to these signalling pathways. One could use KN-93 and autocamtide-2 
related inhibitory peptide [86] for a study to delineate these effects. Our coupled electro- 
physiological and electro-mechanical model also sheds light on the rate dependence of the 
cardiac contractile mechanism. In particular, our model accounts for cAMP-dependent 
modulation of the rate kinetics governing cross -bridge formation. In agreement with 
Janssen [87], we also demonstrate a key linear relationship between the rate of contrac- 
tion and relaxation, which is shown here to be intrinsically coupled over the full range 
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of physiological frequencies both in the absence/presence of ^-adrenergic stimulation. 
This study provides mechanistic, biophysically based explanations for the rate-dependent 
Ca?^ signalling underlying the force-frequency response in rat ventricular myocytes, 
generating useful and testable hypotheses. 

Appendix 

Below is the set of equations modified in the model. 

Equations for currents modified (from Krishna et al. [15]) in the model 
L-Type Ca^"*" current 

Cci}^ current through the DHP- sensitive lca,L channel 

^Ca,L = J^Ca,L ^ ^cAMPJcal (1) 

FcAMPjcal = 1-094 - 0.163 X exp(-0,219x[cAMP] ) (2) 
= 550 + 6.0 X CaMKIIact + CaNact (3) 
^l^ = SSQ + CaMKIIact 

X ( 1.348 + C^™4f ^\ 

V 3.135 X 106 _ 0.755 x CaMKlfiJ^ ) 

§- J^Ca,L described in Krishna et al. [15]. 

t- § described in Krishna et al. [15] causes rate-dependent modulation of the Ica,L 
channel. 

Uptake of Ca^+ from the cytosol into the LSR 

Differential equation for phospholamban phosphorylation 

= kf^B(l + (CaNact X lO-^f)PLBp - ^ff^(l + {CaMKIIact x 10^)^ 
at \p j 

+ ^cAMP,SERCA^^^^dp 

FcAMP,SERCA = 0.1094 - 0.0163 x exp (-0. 219 x[cAMP] ) (6) 

PLBp = l-PLBdp (7) 
/^f^=6800 s-i ; ^1^^=1000 s-^; (Krishna et al. [15]). 



Equations governing electro-mechanics modified (from Rice et al. [55]) in the model 

Regulatory Ca"^"*" -binding to troponin 



dCaTropH , ^ _ . ^ 2+1 x 

— =konTTnIu[ Ca ^\myo (1 " CaTropn) 

dt ^ (8) 

- KffmCaTropH 



dCaTropL =j^^^^Tyii^[Ca^+]myo (1 - CaTropi) 

dt (9) 

- koffirCaTropL 

konT = 22.22 fjLM~^s~^ ; koffHT = 17.36 s"^; koffiT = 173.61 s~^; (Rice et al. [55]). 
dTnIp 

= konTi^PKATnlu - koffTiTnIp (10) 
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KnTi = 698.69 s-^ ; koffn = 80.0 s~^; (estimated from Roof et al. [88]). 

Tnlu = 1 - Tnip (11) 

03x[cAMP] , , 

ApKA = ^ ^ (12) 

^^^^ [c:AMP]+12.1 ^ ^ 

.Force 

(13) 

hfTnew = hfT^ X FcAMP,Force (14) 
,Force (15) 
gxbT„ew = gxbT^ X f cAMP .Force 

(16) 

FcAMP,Force = 1-873 - 1.4 X e;v:/7(-0.192x [cAMP] ) (17) 
Rate constants governing cross-bridge kinetics described in Rice et al. [55]. 
Abbreviations 

[Cc/^+]: calcium ion concentration; [Cc/^+]ys/?: luminal Cc/^+ concentration in the jSR; [Ca^^]myo'- myoplasmic Co^+ 
concentration; [Cc/^+]o: extracellular Cc/^+ concentration; [Ca^^]ryr: Ca^^ concentration at the "mouth" of the RyR 
channel on the dyadic side; CaM: calmodulin; CaMKII: Cc/^+/calmodulin-dependent protein kinase II; CaMKIIact- activated 
Cc/^+/calmodulin-dependent protein kinase II; cAMP: cyclic adenosine monophosphate; CaN: calcineurin; CaNact- 
activated calcineurin; CDF: calcium-dependent facilitation; CDI: calcium-dependent inactivation; CICR: calcium-induced 
calcium-release; CRCP: contraction-relaxation coupling point; DCU: dyadic coupling unit; DHP: dihydropyridine; DHPR: 
dihydropyridine receptor; E-C: excitation contraction; EC^q^: affinity of backward Cc/^+ flux from LSR to myoplasm; EC^^^ = 
affinity of forward Cc/^+ flux from myoplasm to LSR; lca,L- L-type Co^+ current; Icyt.serca- Co^^ uptake current directed from 
the myoplasm to the SERCA; lNa,b'- background sodium current; iNaCa- sodium calcium exchanger current; InoCs- sodium 
cesium pump current; Irmca- plasma membrane Cc/^+ ATPase pump current; Iryr: Ca^^ current due to CICR from an 
individual jSR; lserca,sr'- Ca^^ uptake current directed from the SERCA to the LSR; jSR: junctional portion of the sarcoplasmic 
reticulum; LCC: L-type DHP-sensitive Cc/^+ channel; L-type: long lasting type; LSR: longitudinal portion of the 
sarcoplasmic reticulum; mM: milli molar; mV: milli volt; [Na'^]myo'- myoplasmic A/6I+ concentration; NCX: No'^/Co^^ 
exchanger; nM: nano molar; NO: Nitric oxide; pC: pico coulomb; Pq: Open probability; PKA: protein kinase A; PLB: 
phospholamban; PLB^p unphosphorylated phospholamban; PLBp. phosphorylated phospholamban; PSR: 
phospholamban to SERCA ratio; Ry: ryanodine; RyR: ryanodine receptor; SERCA: sarcoplasmic reticulum Ca^+ ATPase; SL: 
sarcolemma; SR: sarcoplasmic reticulum; VC: voltage clamp; VDI: voltage-dependent inactivation. 
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